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ABSTRACT 

Cosmological N-body simulations were performed to study the evolution 
of the phase-space density Q = p/a 3 of dark matter halos. No significant 
differences in the scale relations Q oc a~ 21 or Q oc M~°- 82 are seen for “cold” 
or “warm” dark matter models. The follow up of individual halos from z = 10 
up to the present time indicate the existence of two main evolutionary phases: 
an early and fast one (10 > z > 6.5), in which Q decreases on the average by a 
factor of 40 as a consequence of the randomization of bulk motions and a late 
and long one (6.5 > z > 0), in which Q decreases by a factor of 20 because of 
mixing induced by merger events. The study of these halos has also evidenced 
that rapid and positive variations of the velocity dispersion, induced by merger 
episods, are related to a fast decrease of the phase density Q. 

Key words: dark matter halos, merger and accretion, phase-space evolution, 
mixing 


1 INTRODUCTION 

Presently, the large-scale structure in the galaxy dis¬ 
tribution on scales 0.02 < k < 0.15h Mpc -1 is success¬ 
fully explained in the cold dark matter (CDM) paradigm 
(Cole et al. 2005 and references therein). The CDM 
power spectrum on these scales derived from large red- 
shift surveys as, for instance, the Anglo-Australian 2- 
degree Field Galaxy Redshift Survey (2dFGRS), is also 
consistent with the Lyman-a forest data in the redshift 
range 2< z <4 (Croft et al. 2002; Viel et al. 2003; Viel, 
Haehnelt & Springel 2004). 

In spite of these impressive successes, there are still 
discrepancies between simulations and observations at 
scales < 1 Mpc, which have been extensively discussed 
in the literature in the past years. In particular, the 
sharp central density cusp in dark matter halos, pre¬ 
dicted by simulations and not seen in the rotation curve 
of bright spiral galaxies (Palunas & Williams 2000; de 
Blok et al. 2001). We could also mention the large num¬ 
ber of subhalos present in simulations but not observed 
(Kauffmann, White & Guiderdoni 1993; Moore et al. 
1999; Klypin et al. 1999), as in the case of our Galaxy 
or M31. Moreover, deep surveys (z> 1-2) as the Las 
Campanas Infrared Survey, HST Deep Field North and 
Gemini Deep Deep Survey (GDDS) are revealing an ex¬ 


cess of massive galaxies with respect to predictions of 
the hierarchical scenario (Glazebrook et al. 2004). 

Massive halos are formed by continuous accretion 
and/or by merging. After a merger episode, the result¬ 
ing halo is not in equilibrium, but after few dynamical 
time scales (tdyn ~ (Gp) _1,/2 ) the gravitational poten¬ 
tial becomes almost steady. Violent relaxation produces 
a more mixed system (Tremaine, Henon & Lynden-Bcll 
1986), reducing the value of the coarse-grained distri¬ 
bution function (DF), for which Q (defined as the ra¬ 
tio between the density and the cube of the 1-D veloc¬ 
ity dispersion in a given volume, e.g., Q = p/a 3 ) is an 
estimator. Analyses of the core phase-space density of 
dwarf spheroidal galaxies, rotating dwarfs, Low Surface 
Brightness (LSB) galaxies and clusters of galaxies sug¬ 
gest that Q oc a~ n , with n ~ 3-4 (Dalcanton & Hogan 
2001). This behavior can be understood if the merging 
halos were initially almost in equilibrium and the fusion 
process preserves approximately the physical density as 
each layer is homologously added to form the new sys¬ 
tem (Dalcanton & Hogan 2001). A scaling relation close 
to Q oc a~ 3 was obtained from cosmological simula¬ 
tions by Dave et al. (2001) for collisionless dark matter 
as well as for self-interacting dark matter (SIDM). The 
result for SIDM is unexpected since in this case the ma¬ 
terial should be compressed to higher densities during 
the merger event, sinking to where it reaches the lo- 
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cal pressure equilibrium and where the specific entropy 
matches. Therefore the merger should occur at nearly 
constant Q rather than constant density. Besides such a 
scaling, high resolution simulations of galaxy-size CDM 
halos indicate an increase of the coarse-grained phase- 
space density Q towards the center (Taylor & Navarro 
2001). Similar results were obtained by Rasia, Tormen 
& Moscardini (2004), who have also obtained a power 
law variation, e.g., Q tx r _/3 for cluster-size halos, with 
(3 quite close to the value found by Taylor & Navarro 
(2001), namely, f3 « 1.87. These results suggest that 
the density profiles of dark matter halos are a conse¬ 
quence of a hierarchical assembly process that preserves 
the phase-space stratification (Taylor & Navarro 2001) 
or even a generic feature of violent relaxation (Williams 
et al. 2004). 

In the present paper we report new cosmological 
simulations aiming to study the evolution of the “poor 
man’s” phase-space density Q in the core of dark mat¬ 
ter halos. In particular, we focus on the role of accre¬ 
tion/merger processes, which affect the dynamical re¬ 
laxation of halos, leading to structures with higher en¬ 
tropy or, equivalently, with lower phase-density Q. We 
show that the Q-cr scaling independs practically on the 
mechanism by which halos acquire mass. Notwithstand¬ 
ing, a small differencial effect seems to exist, since halos 
which have undergone major merger events have slightly 
lower phase-space densities than halos of same mass, 
which have accreted matter quietly. The same behav¬ 
ior is noticed if warm dark matter is considered instead 
of CDM. The second part of this investigation concerns 
the follow up of some individual halos in order to obtain 
details of changes in the phase-space density, occurring 
during merging events. We have found an early rapid de¬ 
crease of the phase-density Q, during and just after the 
first shell crossing, followed by a slow and long decrease 
as halos accrete mass. Particles of captured subhalos 
diffuse in the velocity space and after few dynamical 
timescales, share with particles of the main halo a com¬ 
mon (Gaussian) velocity distribution. This paper is or¬ 
ganized as follows: in Section 2 the N-body simulations 
are brifly described, in Section 3 the evolution of the 
phase-space density Q as well as the effects of merger 
events are discussed and, finally, in Section 4 the main 
conclusions are given. 


2 NUMERICAL SIMULATIONS 

In this work, we use the same halo catalogs as in 
Peirani, Mohayaee & de Freitas Pacheco (2004, here¬ 
after PMP04), but in the present study objects with 
masses higher than 10 13 Mq were also included. For the 
sake of completeness, we summarize here the main steps 
performed to prepare these catalogs. 

The N-body simulation uses the adaptive particle- 
particle/particle-mesli (AP 3 M) code HYDRA (Couch- 
man, Thomas & Pearce 1995). The adopted cosmo¬ 
logical parameters were li = 0.65, = 0.3 and 17 a 

= 0.7, with the power spectrum normalization erg = 
0.9. The simulation was performed in a box of side 
301V 1 Mpc including 256 3 particles, corresponding to 


a mass resolution of 2.05 x 10 8 Mq. The simulation 
started at z = 50 and ended at the present time (z = 
0). A similar simulation (here referred as S2) was also 
performed for a “warm” dark matter (WDM) model. 
In this case, the fluctuation power spectrum becomes 
P(k) = AkT^ drn (k)T^ drn (k), where the adopted transfer 
function for WDM is given by (Bode, Ostriker & Turok 
2001) 

T wd m(k) = [1 + (afe) 2 ’ 4 ] _25/6 (1) 

where a = 0.102 for WDM particles having a mass of 
0.5 keV and k is in hMpc -1 . This particle mass corre¬ 
sponds to a free-streaming mass of about 1.7 x 10 11 Mq, 
scale below which structures become to be suppressed. 
A third simulation (here referred as S3) with the same 
cosmological parameters but including 128 3 CDM parti¬ 
cles was also run. However, particles were separated into 
two categories, the first having masses ten times higher 
than the second. This simulation was carried out using 
the code GADGET (Springel, Yoshida & White 2001) 
and was used to study the energy transfer between par¬ 
ticles during merger events. Here, the “softening” length 
was set equal to 12 h _1 kpc. 

Halos were initially detected by using a friends- 
of-friends (FOF) algorithm and, in a second step, un¬ 
bound particles were removed by an iterative proce¬ 
dure. Thus, halos of both catalogs are all gravitationally 
bound objects. For further details, the reader is referred 
to PMP04. 

The accretion catalog, comprising 781 objects, in¬ 
cludes halos which have never undergone a major 
merger event and whose masses varied continuously and 
smoothly. The merger catalog contains 567 halos which 
had at least one major merger episode in their history, 
corresponding to an increase of their masses at least by 
a factor of 1/3 in the event. In latter class, violent vari¬ 
ations of the gravitational potential occur while in the 
former category the potential varies smmothly all the 
time. Halos were followed from z = 3.5 until z = 0, but 
about 40 objects, with enough particles (N > 50), were 
followed during a longer time interval (from z = 10 until 
z = 0). From the simulation S2, 216 halos were included 
in the merger catalog and 334 halos in the accretion cat¬ 
alog (3.5> 2 > 0) but only 5 halos have been followed 
in the interval 10 > 2 > 0. 


3 THE PHASE-SPACE DENSITY 
EVOLUTION 

3.1 Scaling relations 

One of the goals of this work is to study the evolution 
of the phase-space density Q = p/o 3 , in order to verify 
how the “violence” or the “quietness” of the accretion 
process may affect the relaxation road to equilibrium 
and the resulting scaling relations Q tx a~ n or Q oc 
M~ k . 

For each halo at a given redshift, the value of Q was 
computed by adopting the following procedure. Firstly, 
the halo center is defined as being coincident with the 
position of the most bound particle, since during a 
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Figure 1. The phase-space density Q for halos of both cata¬ 
logs as a function of the ID velocity dispersion a is shown for 
the CDM model (upper panel) and for WDM model (lower 
panel). In each panel, the solid line indicates the best fit solu¬ 
tion (Q oc a 2 2 ). For comparison, lines representing Q oc < 7 " 3 
(dotted) and Q oc a ~ 4 (dashed), normalized at the same 
point, were also included in the plot. 


merger event the center of mass is “ill-defined”. In a sec¬ 
ond step, all particles inside a radius r* around the cen¬ 
ter are selected and r* is taken to be equal to one tenth 
of the gravitational radius, e.g., r* = 0.1GM 2 /| W |, 
where W is the total gravitational energy of the system. 
This radius was chosen in order to include enough par¬ 
ticles necessary to compute adequately the density as 
well as the velocity dispersion and still be representa¬ 
tive of the core region. Instead of evaluating the core 
density by a simple mean ( p c = 3mN c /Anr J), we have 
estimated the density pi at the position of each selected 
particle i by using the Voronoi tessellation procedure in 
order to increase the accuracy of the density estimation 
when the number of particles is small. In this technique, 
the 3-D space is divided into polyhedral cells centered 
on each particle, with boundaries defined in order to in¬ 
clude all points closer to the center (see, for instance, 
van de Weygaert 1994, for details). The cell density is 
simply the mass of the particle divided by the volume 
of its cell. Calculations were made with the help of the 
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Figure 2. The phase-density Q as a function of the core mass 
for halos (CDM model) of both catalogs at z = 0. The solid 
line shows the best fit (Q oc m~ 0S3 ) and the small panel 
shows the same plot, but discriminates objects of different 
catalogs. The vertical bar indicates a typical uncertainty in 

Q- 


public package QHULL*. Finally, the core density of the 
halo is estimated from 


Pc = (1/JVc) 


( 2 ) 


and the 1-D velocity dispersion from 
o = \J (ui + o 2 + erf)/3. 

In Fig. 1, the phase-space density Q for all halos 
and for both CDM and WDM models is shown as a 


function of the 1-D velocity dispersion a at z — 3.0 and 
z = 0. For the CDM model, a best fit to the data indi¬ 
cates that Q oc ct~ 2 08 for z = 3.0 and Q oc cr~ 216 for 
z = 0. No significant modifications in the exponent n 
( Q oc o~ n ) are noticed if the definition of r* is changed. 
For instance, at z = 0, n = 2.23 and n = 2.18 if r* is 
decreased by a factor of two or increased by 50% respec¬ 
tively. We have also investigated how our results could 
be affected by using the mean density p c instead of the 
tessellation method. In fact, a slight increase in the ex¬ 
ponent n is observed (n = 2.64 at z = 3 and n = 2.53 
at z = 0). The reason for such differences is the follow¬ 
ing: the density evaluated from the Voronoi procedure 
in a given volume is weighted by the inverse cell volume 
associated to each considered particle. Thus, if impor¬ 
tant density gradients are present, as those observed 
in DM halos, the two density estimations differ. When 
using simulated data, this effect in more pronounced 
when massive halos are considered, since the number 
of particles is large, strengthening the gradient effect. 
As a consequence, the use of a simple mean underesti¬ 
mates the density as well as the phase-density Q with 
respect to the Voronoi method. This effect is minimized 
for low mass halos, since the number of particles inside 


* http://www.geom.umn.edu/software/qhull/ 
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the considered volume is small, diluting the gradient ef¬ 
fect. In order to give some numbers, the core volume 
of a halo of 7 x 10 12 M© including ~ 3200 particles is 
comparable to the total cell volume, but the density 
derived by the tessellation procedure is about a factor 
of two higher than the simple mean density. However, 
both density estimations give similar results if a typi¬ 
cal halo of 1.2 x 10 11 M©, including ~ 55 particles, is 
considered. Some differences can also be seen when our 
simulated data is compared with observations (see, for 
instance, Fig. 1 in Dalcanton & Hogan 2001), since Q 
values estimated for low brightness galaxies are higher 
than our expectations. This is probably a consequence 
that density estimations from observational data include 
the contribution of baryonic matter, dominant in the 
central regions of these objects (Ortega & de Freitas 
Pacheco 1993), which is not the case for clusters. 

It is worth mentioning that the exponent n found 
from our simulations can be understood on the basis of 
the following considerations: from data of our catalogs, 
halos acquire mass continuously and evolve according 
to the scaling relation M oc a s ^ 3 . If they evolve almost 
in virial equilibrium in the considered redshift interval, 
then Q oc in good agreement with our results. It 

should also be emphasized that by the time of halo for¬ 
mation, the particle thermal velocities are quite small 
and the early random motions, present when simula¬ 
tions begin, are a consequence of gravitational “heating” 
by primordial density fluctuations. This is true whether 
cold or warm dark matter is considered. Thus, if scaling 
relations can be checked by simulations, this is not the 
case for the normalization of the initial phase-space den¬ 
sity. Moreover, as halos evolve from z = 3 to z = 0, the 
Q-<t relation becomes slightly but significantly steeper 
and, at a given velocity dispersion, the value of Q de¬ 
creases by a factor of 4 in the same redshift interval. 
A recent perturbative analysis of the spherical gravita¬ 
tional collapse (Amarzguioui & Grocjtn 2005) indicates 
that the entropy increases (e.g., the phase-space Q de¬ 
creases) in the process, as a consequence of potential 
energy conversion into random kinetic motions. This is 
in agreement with our findings, which covers the non¬ 
linear regime. 

Halos in a WDM model have, on the average, sys¬ 
tematically higher virial ratios (2 T/ \ W |) than CDM 
halos, i.e., they are less relaxed. In spite of this, no signif¬ 
icant differences are seen in the Q — a diagram between 
both dark matter models. The exponent is equal to n = 
1.96 (z = 3.0) and n = 2.13 (z = 0), comparable to the 
values found for the CDM model. 

The evolution of Q as a function of the core mass of 
the halo for objects included in the accretion and in the 
merger catalogs is shown in Fig. 2, for 2 = 0 in the case 
of CDM (similar results are obtained for WDM). Data 
were billed in logarithm intervals of A log(A4 /Mq) = 
0.20 for a higher clarity of the plot. Best fit solutions 
give Q oc M ~ 0 80 for the accretion and Q oc A/ -0 ' 82 for 
the merger samples respectively. Again, using the em¬ 
pirical scaling relation derived from our simulated data, 
e.g., M oc a 8 / 3 and the assumption that halos evolve 
in quasi-equilibrium from 2 = 3.0 up today, one obtains 


M oc a 7 ^ 8 , quite comparable to the values derived from 
simulated data. 


3.2 Study of individual halos 

Halos are continuously accreting mass (violently or not). 
This process affects their dynamical state and con¬ 
tributes to decrease their phase-space density Q. Thus, 
scaling relations as Q oc o~ n or Q oc M~ k should be 
considered in a statistical sense and they are not ex¬ 
pected to describe the evolution of any particular halo, 
whose evolution depends on his accretion history. 

In order to understand the early behavior of the 
coarse-grained DF (measured by Q), 43 CDM and 5 
WDM halos in the mass range M ~ 10 11 ” 12 M© were 
followed since z = 10 up today. Criteria to identify a 
given halo at different redshifts were discussed in our 
previous work (PMP04). For each halo, we have evalu¬ 
ated the phase-space density Q, the core density p c and 
the velocity dispersion a as a function of time. Some 
typical examples are exhibited in Fig. 3, where the evo¬ 
lution of these parameters for 3 individual halos (2 CDM 
and one WDM) is shown. 

As halos accrete mass, they approach a state of 
dynamical equilibrium, measured by their virial ratio 
2T/| W |, which approaches asymptotically the unity. 
On the one hand, their gravitational radius increases 
continuously and, as a consequence, the core density p c 
decreases on the average by one order of magnitude in 
the interval 10< 2 < 0, for halos in the considered mass 
range, indicating that the evolution of Q can not be de¬ 
rived by a simple homologous model. On the other hand, 
from z = 10 up to about z ~ 6.5, there is a rapid increase 
of the velocity dispersion, associated to the process of 
relaxation, followed by a phase of slow “heating”. In the 
early evolutionary phases, an important energy transfer 
from bulk to random motions occurs due to collective 
effects, heating all the particles, regardless their initial 
energies. Such a transfer by collective effects is clearly 
seen in merger episodes issued from simulation S3, in¬ 
dicating that particles with different masses have the 
same final dispersion velocity. This relaxation process 
reduces considerably the phase-space density in a short 
time interval. Q decreases, on the average, by a factor 
of 40 in a time interval of about 0.5 Gyr, whereas in the 
next phase, which lasts about 13 Gyr, the phase-space 
density decreases by factor of 20 due to phase-mixing 
mechanisms induced by mass accretion processes. 

The effect of the “heating” in variations of the 
phase-space density can be clearly seen if the time 
derivatives of Q and a are compared. The time deriva¬ 
tive of the velocity dispersion has, in general, clear peaks 
indicating phases of fast “heating”. The first and princi¬ 
pal peak is associated to the “randomization” of initial 
bulk motions and the others are associated to merger 
events. The amplitude of these secondary peaks depends 
on the mass of the captured subhalo. These events pro¬ 
duce a sudden decrease of the coarse-grained DF and 
a more relaxed system, with a higher entropy. In Fig. 
4 we show three examples of CDM halos where effects 
of merger events seen in the evolution of the velocity 
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Figure 3. The evolutions of Q, a and p c for 3 halos from z = 10 to z = 0. Panels in the first two columns concern halos 
evolving in a CDM model while panels in the last column concerns an example of a halo evolving in a WDM model. 


dispersion are clearly anti-correlated with phase-space 
density variations. 


4 CONCLUSIONS 

The dynamical evolution of dark matter halos depends 
on their accretion history. During merger events strong 
mixing effects occur away from an equilibrium state, 
driven by large-scale fluctuations of the gravitational 
potential. In this work we have investigated different 
aspects of mixing through analyses of the phase-space 
density Q evolution after accretion/merger events. 

Our simulations indicate that for 2 < 3.0 the phase- 
space density scales as Q tx a~ 21 and that such a scal¬ 
ing does not depend on the adopted dark matter model. 
The same behavior is obtained either for cold or warm 
dark matter, in spite of CDM halos be always more re¬ 
laxed, at a given redshift, than WDM halos of com¬ 
parable mass. From the scaling relation derived from 
our simulations between the halo mass and the velocity 
dispersion, considering also that at those redsliifts halos 
have virial ratios near unity, we obtain Q oc a'^ 3 , consis¬ 
tent with the results derived directly from our numerical 
experiments. 

Q values derived for halos which have underwent 
important merger events are, on the average, slightly 


lower than those derived for halos which have grown 
by accretion only. This is consistent with the expecta¬ 
tion that violent variations of the gravitational poten¬ 
tial, which occur in a merger event, lead to a more mixed 
system. Additional information about the evolution of 
the phase-space density Q was obtained from the follow 
up of individual halos from 2 « 10 up today. In the early 
evolutionary phases (10 > 2 > 6.5), Q decreases on the 
average by a factor of 40, as a consequence mainly of 
the randomization process of the initial bulk motions. 
In a late (and long) phase (6.5 > z > 0), Q decreases by 
a factor of 20 on the average, consequence of continu¬ 
ous mass accretion. These processes are well evidenced 
when the time derivative of the phase density Q and that 
of the velocity dispersion for individual halos are com¬ 
pared. There is a clear anticorrelation between these two 
quantities: rapid positive variations of the velocity dis¬ 
persion, induced by merger events, are connected with 
a fast decrease of the phase density Q. 
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Figure 4. The evolution of the time derivative of Q and those of the velocity dispersion for 3 examples of CDM halos. The first 
event labelled 1 corresponds in general to the randomization of bulk motions, while the others correspond to merger events. 
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